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ABSTRACT 

We present results from a blind survey to identify strong gravitational lenses among the pop¬ 
ulation of low-redshift early-type galaxies. The SINFONI Nearby Elliptical Lens Locator 
Survey (SNELLS) uses integral-field infrared spectroscopy to search for lensed emission line 
sources behind massive lens candidates at z< 0.055. Lrom 27 galaxies observed, we have 
recovered one previously-known lens (ES0325-G004) at z =0.034, and discovered two new 
systems, at z = 0.031 and z = 0.052. All three lens galaxies have high velocity dispersions 
(tr > 300kms“^) and a-element abundances ([Mg/Le] >0.3). Lrom the lensing configura¬ 
tions we derive total J-band mass-to-light ratios of 1.8±0.1, 2.1±0.1 and 1.9±0.2 within the 
~2 kpc Einstein radius. Correcting for estimated dark-matter contributions, and comparing to 
stellar population models with a Milky Way (Kroupa) initial mass function (IML), we deter¬ 
mine the “mass excess factor” , a. Assuming the lens galaxies have “old” stellar populations 
(10±1 Gyr), the average IML mass factor is (a) = 1.10±0.08±0.10, where the first error is 
random and the second is systematic. If we instead fit the stellar populations from 6dL op¬ 
tical survey spectra, all three galaxies are consistent with being old, but the age errors are 
3-4- Gyr, due to limited signal-to-noise. The IML constraints are therefore looser in this case, 
with (a) = 1.23^0 13 ± 0.10. Our results are thus consistent with a Kroupa IML (a= 1.00) on 
average, and strongly reject very heavy IMLs with a > 2. A Salpeter IML (a = 1.55) is incon¬ 
sistent at the 3.5(7 level if the galaxies are old, but cannot be excluded using age constraints 
derived from the currently-available optical spectra. 

Key words: gravitational lensing: strong — galaxies: elliptical and lenticular, cD — galaxies: 
stellar content — stars: luminosity function, mass function 


1 INTRODUCTION 

Among all techniques for determining galaxy masses, strong gravi¬ 
tational lensing yields the most precise and accurate measurements, 
with errors of only a few per cent in favourable cases (Treu 2010; 
Courteau et al. 2014). With no dependence on the dynamical or 
thermal state of tracer material, lensing is free from many of the 
degeneracies afflicting other mass probes. Being sensitive to all 
gravitating matter, lensing provides some of the strongest direct 
evidence for dark matter halos around elliptical galaxies. If the 
contributions of dark and baryonic mass can be distinguished, e.g. 
through their different spatial distribution, then lensing can be used 
to measure the stellar mass-to-light ratio T. Such measurements 
can provide stringent constraints on stellar population models, and 


* Based on observations collected at the European Southern Observatory, 
Chile ESO Programmes 093.B-0193 and 293.B-5026. 
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especially on the initial mass function (IMF), since faint low-mass 
stars (and/or dark remnants of high-mass stars) add substantially to 
the total mass. 

Analysing the lensing properties and stellar velocity disper¬ 
sions ((t) of a sample of 56 strong-lens systems from the SLAC^ 
survey, Treu et al. (2010) found that the IMF and the dark-matter 
halo profile could not both be universal. Holding the halo pro¬ 
file fixed, they derived a trend in “excess” stellar mass, which 
they attributed to a non-universal IMF. Galaxies with lower veloc¬ 
ity dispersion were consistent with Milky-Way-like (e.g. Kroupa 
2001 or Chabrier 2003) IMFs, while for their most massive galax¬ 
ies, with (T > 300kms“^, the SLAGS modelling requires stellar 
mass-to-light ratios twice that expected for a MW-like IMF, i.e. 
Q = T/Tref ~2 (where Tref is the expected stellar mass to light 
ratio for a Kroupa IMF, given observational constraints on the star- 


^ Sloan Lenses Advanced Camera for Sui-veys. 
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formation history). Auger et al. (2010) later showed that heavier- 
than-MW IMFs were still required on average even when allowing 
for contraction of halos in response to the stellar component, fol¬ 
lowing Gnedin et al. (2004). With this modification however, the 
requirement for a trend in a with galaxy mass is not significant, 
and even the most massive lenses have a « 1.4, intermediate be¬ 
tween Kroupa (2001) (a = 1) and a Salpeter IMf0 (a = 1.55). The 
difference with respect to the Treu et al. result highlights the chal¬ 
lenge in separating baryonic from dark-matter contributions to the 
total lensing mass. 

Concurrently with the SLACS analysis, two other largely in¬ 
dependent methods also found evidence to favour heavy IMFs in 
giant ellipticals. Fitting dynamical models to integral-field stel¬ 
lar kinematics from the SAURON survey, Cappellari et al. (2006) 
showed that excess mass was required in the most massive ellipti¬ 
cals, which could be attributed either to an increasing dark matter 
fraction or to a systematic trend in the IMF. With more detailed 
modelling of the larger ATLAS3D survey dataset, Cappellari et 
al. (2012, 2013) were able to reject the dark-matter solution, for a 
range of possible halo prescriptions (see also Thomas et al. 2011). 
They concluded that the IMF mass normalisation indeed increases 
with velocity dispersion, though the derived trend has a shallower 
slope than that of Treu et al. (2010). At 300kms“^ the (extrap¬ 
olated) ATLAS3D trend implies an average a= 1.6, with an esti¬ 
mated intrinsic scatter of 20 per cent. The two trends are not nec¬ 
essarily inconsistent however, since the ATLAS3D sample probes 
smaller velocity dispersions than the SLACS lenses; Posacki et al. 
(2015) derived a quadratic a-vs-a relation, with the slope steepen¬ 
ing at high (7, which adequately describes both samples. An alter¬ 
native approach to dynamical modelling using single-fibre spectra 
for large galaxy samples similarly favours an IMF trend reaching 
a > 2 at a > 300kms“^ (e.g. Dutton et al. 2013; Conroy et al. 
2013; Tortora et al. 2014). 

In parallel with these works, van Dokkum & Conroy (2010) 
studied gravity-sensitive spectral features in giant ellipticals, and 
found evidence for an excess population of cool dwarf stars, cor¬ 
responding to a bottom-heavy IMF. The spectroscopic method is 
quite distinct from lensing or dynamics, in that the stellar content 
is probed directly, rather than through any constraint on the total 
mass, and is hence independent of uncertainties due to dark matter. 
Further studies with this approach (Spiniello et al. 2012; Conroy 
& van Dokkum 2012b; Smith, Lucey & Carter 2012; Ferreras et 
al. 2013; La Barbera et al. 2013; Spiniello et al. 2014) confirmed 
a trend of increasing dwarf-star content as a function of galaxy 
“mass” (velocity dispersion or correlated quantities, e.g. element 
abundances). Conroy & van Dokkum (2012b, CvD12b) analysed a 
sample of ~30 galaxies and derived an implied mass excess (as¬ 
suming a three-part power-law IMF) of a«2 at CT = 300kms“'^, 
similar to the SLACS results. Recent improvements of the models 
and enlargements to the sample, however, yield a wider galaxy-to- 
galaxy range in a at high mass, and hence a smaller average a, but a 
Milky-Way-like IMF remains disfavoured (Conroy & van Dokkum, 
in preparation). 

This paper describes an effort to derive robust constraints on 
a in high velocity dispersion galaxies, by finding and analysing 
strong-lensing galaxies at much lower redshift than those in exist- 


^ As usual in extragalactic studies, we (mis)use the term “Salpeter IMF” to 
refer to a single-slope power law At(M) oc with a slope of x = 2.35 

over the wide mass range 0.1-100 Mq. Salpeter (1955) derived this slope 
only over the range 0.5-10 Mq . 


ing samples. The accurately determined quantity in a strong-lens 
system is usually the total projected mass within the Einstein aper¬ 
ture, defined by a threshold in enclosed surface density (which de¬ 
pends on the angular diameter distances of the lens and source). 
For a given source redshift, the Einstein radius i?Ein (in angular 
units) varies only slowly with lens redshift. Hence, while distant 
lens galaxies (e.g. the “Cosmic Horseshoe” studied by Spiniello et 
al. 2011) have i?Bin comparable to, or larger than, the half-light ra¬ 
dius (Res) of the stellar component, in nearby lenses the mass can 
be determined at a small fraction of ReS- Since stars are more cen¬ 
trally concentrated than dark matter, nearby lenses are optimal for 
measuring the stellar mass-to-light ratio and hence for constraining 
the IME in massive galaxies. A further motivation for studying such 
systems is that the lens galaxies can be characterised in great detail, 
with high signal-to-noise spectroscopy, to determine the probable 
star-formation history, which in turn tightens the limits on a. 

Only one low-redshift (here defined as 2 : <0.05) massive 
elliptical strong-lens galaxy has been reported before novj^. 
ES0325-G004 is a cr = 330km s“^ galaxy at 2 = 0.034 which 
lenses a 2 = 2.14 background source into a partial Einstein ring with 
F?Bin = 2.85 arcsec (Smith et al. 2005; Smith & Lucey 2013, here¬ 
after SL13). Analysing the ES0325-G004 system, SL13 derived a 
stellar mass-to-light ratio fully consistent with a Milky Way IME, in 
apparent contradiction to the evidence for heavy IMEs in massive 
galaxies from SLACS, CvD12b and ATLAS3D. 

A key question is whether lensing masses for giant ellipticals 
in general favour a MW-like IME or if ES0325-G004 is simply 
an “outlier” in a broad distribution in a. To construct a mean¬ 
ingful sample of low-redshift strong-lensing early-type galaxies, 
as required to address this issue, we are conducting the SINEONI 
Nearby Elliptical Lens Locator Survey (SNELLS). Briefly, our ap¬ 
proach is to search within infra-red integral-field datacubes, to de¬ 
tect background line-emitters behind a carefully-selected sample of 
candidate lenses. In this paper, we give a broad description of the 
survey and present the first results from SNELLS, including the 
identification of two new multiple-image 2 < 0.05 lens systems, 
and limits on the IME in these, plus ES0325-G004, using a uni¬ 
form dataset and methodology. Sectionj^discusses the motivation 
for the survey strategy and describes its practical implementation 
so far. Section[3presents the newly-discovered lenses, and the re¬ 
covery of the ES0325-G004 system. In Section|4] we make a first 
estimate for the lensing masses, mass to light ratios, and the IME 
a factor. The results are compared to previous work in Section 0 
while the conclusions and future prospects are summarized in Sec¬ 
tion 

Eor consistency with other recent work in the field, we 
adopt cosmological parameters from WMAP7 where necessary, i.e. 
F/o = 70.4kms“^ Mpc“^, Dm = 0.272 and Da =0.728 (Komatsu 
et al. 2011). We comment in Section l43] on the effect of using in¬ 
stead the Planck Collaboration (2014) cosmological parameters. 


® The lens of the Einstein Cross system at 2 = 0.04 is a spiral galaxy. The 
possible extended arc around in NGC5532 (Hardcastle et al. 2005) has not 
yet been confirmed with deeper imaging or spectroscopy. Several cluster- 
scale, rather than galaxy-scale, arc candidates have been reported at low 
redshift (Campusano, Kneib & Hardy 1998; Blakeslee et al. 2001). The 
nearest objects in SLACS are at 2 ~ 0.06; they are among the least massive 
in the sample and do not require heavy IMFs. 
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2 A BLIND SURVEY FOR LENSES AMONG 
LOW-REDSHIFT ELLIPTICALS 

The two principal obstacles to assembling a sample of low-redshift 
strong-lensing massive ellipticals are the intrinsic rarity of the lens 
population and the low contrast between lens and source in “typi¬ 
cal” systems. Any search restricted to very bright arcs will have a 
low yield of lenses per candidate, and hence must probe a very large 
survey volume (e.g. Gpc® within 2 ; < 0.3 for SLACS). To find 
lenses in useful numbers within a “local” volume (say 0.07 Gpc® 
within 2 < 0.06), a search must be able to detect systems in which 
the lensed source is very faint, a task made more difficult by super¬ 
position on a bright foreground galaxy. Traditional methods strug¬ 
gle to identify lenses in this regime. 

As we show in this section, targetted integral field unit (IFU) 
observations, especially in the infra-red (IR), can provide the im¬ 
proved contrast necessary to discover less visually-spectacular, but 
more numerous, lens systems at low redshift. 

2.1 IFU lens-search strategy 

Previous optical searches for galaxy-scale strong lenses fall mainly 
into two categories, each with its strengths and limitations: 

Spectroscopic approaches, like SLACS (Bolton et al. 2006) 
look for anomalous emission lines from background galaxies in in¬ 
tegrated spectra from large surveys (see also Warren et al. 1996; 
Hall et al. 2000). Because the spectrum is obtained from a large 
aperture (a 3” fibre for the Sloan Digital Sky Survey, SDSS), only 
very bright-lined sources can be detected over the dominant con¬ 
tinuum flux from the lens galaxy. In SLACS, moreover, the SDSS 
spectral range imposes a limit on the source redshift (2src ^ 1). 
and the search volume behind each target is correspondingly small. 
After spectroscopic selection, follow-up data is needed to establish 
the image positions, either using HST to heighten contrast against 
the lens galaxy, or with an IFU to map the line emission. 

Morphological methods (e.g. Cabanac et al. 2007) instead 
search systematically for arc-like features in broad-band images, 
either using automated algorithms (e.g. Belokurov et al. 2009; 
Gavazzi et al. 2014) or visual inspection (including crowd-sourced 
classification, e.g. Marshall, Lintott & Fletcher 2014). When lensed 
images are found, follow-up spectroscopy is necessary to estab¬ 
lish the source and lens redshifts, before physical properties can be 
determined. With ground-based data, this approach becomes less 
effective for detecting low -2 lenses, where the arcs lie inside the 
effective radius of the lens galaxy and are hard to distinguish from 
stellar shells, spiral arms, tidal features, etc. 

The rarity of high-contrast lenses can be overcome by apply¬ 
ing either the spectroscopic or the morphological search techniques 
to large (and pre-existing) survey datasets, but most of the result¬ 
ing lenses are then inevitably found at moderate redshifts. Neither 
method is well suited to detecting faint arcs behind nearby massive 
galaxies. 

ES0325-G004 was serendipitously identified as a lens 
through inspection of Hubble Space Telescope images (Smith et al. 
2005). In this case, the necessary contrast between lens and source 
is provided by the high angular resolution of the images, and by 
the much bluer colour of the arcs compared to the elliptical. Per¬ 
haps surprisingly, archival HST data is not presently available for 
many very massive nearby ellipticals, except for those at the centres 
of rich clusters, which are poor candidates for galaxy-scale lenses 
(see Section l2.2.1l l. The discovery of ES0325-G004 suggests that 
a systematic blind search with HST for arcs behind candidate lenses 


might be fruitful, but it would be an observationally expensive pro¬ 
gramme, and still require spectroscopic follow up for redshift mea¬ 
surements. 

For the lens search described here, we instead use ground- 
based IFU spectroscopy to combine some of the advantages of both 
the spectroscopic and imaging-based methods. The IFU data con¬ 
centrate the emission spatially as well as in wavelength, increasing 
contrast over the lens galaxy (>40x lower lens galaxy flux in a 
1 " aperture at 0.25xi7eft, compared to a full 3” SDSS fibre). To 
maximize the background volume available per lens candidate, we 
work in the near IR; this also allows us to target the bright Ha 
emission line. A combination of J and H bands probes lensed Ha 
over a volume-per-candidate times larger than that available to 
SLACS (mainly lensed [O III] emitters in optical spectra). Finally, 
since we perform dedicated observations, rather than exploiting an 
existing survey, we maximize efficiency by pre-selecting a sample 
of the most promising lens candidates based on velocity disper¬ 
sion. The lensing cross-section scales as a*‘, for isothermal mass 
profiles, giving a factor of 20 gain at {a) — 315 km compared 
to (crsDSs) = 150kms“^. Together, these factors yield a search 
efficiency which is > 1000 times larger, on a per target basis, than 
the SLACS approach, so that a useful sample of lenses can be as¬ 
sembled from observations of only a few tens of candidates. 

2.2 Implementation 

The specific implementation of our approach uses the SINFONI 
instrument (Eisenhauer et al. 2003; Bonnet et al. 2004), a near- 
IR image-slicing IFU spectrograph at the ESO 8.2m Very Large 
Telescope (VLT), Unit Telescope 4. 

We established the optimal survey design based on the follow¬ 
ing assumptions; 

(i) The lenses have singular isothermal ellipsoid profiles, with 
(jsiE =315 2iens = 0.04 and ellipticity of 0.3. For this 
model, the cross-section for multiple-image lensing is 25 arcsec^ 
for a distant source (2src ^ 0.5). 

(ii) The sources are Ha emitters drawn from the evolving lumi¬ 
nosity function as determined from wide-field narrow-band imag¬ 
ing by Sobral et al. (2013). We adopt a fixed faint-end slope ana = 
— 1.6, a linear relationship for the break luminosity log with 2 , 
and interpolate the evolution of the luminosity density, log 

with a cubic spline. 

(iii) The line detection limit is computed as a function of wave¬ 
length using the ESO exposure time calculator, including the ef¬ 
fects of instrument efficiency, sky emission lines, atmospheric ab¬ 
sorption, etc. The total survey yield is computed by integrating the 
counts down to the detection limit and summing over the whole 
wavelength coverage of the observation. 

In principle, lensing magnification increases the total source flux, 
but initial detection of peaks will not be aided fully by this factor, 
because the signal is also spread into more pixels. Since detectabil¬ 
ity is mainly governed by surface brightness, we conservatively ig¬ 
nored the magnification effect in our calculations. 

Experimenting with the parameters of this calculation shows 
that for a given total exposure time, the yield is maximized by 
splitting the observation between J and H bands, to increase the 
background volume surveyed. Likewise, the slope of the source 
counts function is such that shallow observations of more targets 
are preferable to deep observations of fewer candidates. The final 
predicted survey yield, i.e. identified lenses per observed candidate, 
is ~ 10 per cent, for integration times of -^1 hr split between the J 
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Table 1. Summary of the 27 SNELLS target galaxies observed in ESO Period 93. Redshifts, z, and velocity dispersions, cr are from the 6dE Peculiar Velocity 
Survey (Campbell et al. 2014) or from SDSS DRV (Abazajian et al. 2009) as indicated. 


Target 

NED ID 

source 

Z 

cr [kms 

Lens ID 

Notes 

J0058513-162809 

NGC 0333B 

6dE 

0.0541 

313 ± 20 



J0141423-073528 

2MASX J01414232-0735281 

6dE 

0.0519 

320 ± 18 

SNL-2 


J0145535-065608 

IC0158 

6dE 

0.0523 

328 ± 20 



J0202308-505553 

ESO197-G018 

6dE 

0.0215 

333 ± 24 



J0213108-062026 

MCG-01-06-071 

6dE 

0.0433 

308 ± 21 



J0258177-044906 

MCG-01-08-023 

6dE 

0.0308 

324 ± 27 

SNL-3 


J0302383-603219 

2MASX J03023835-6032198 

6dE 

0.0543 

332 ± 26 



J0420410-144651 

IC 0367 

6dE 

0.0308 

303 ± 20 



J0837138-032141 

CGCG 004-087 

6dE 

0.0402 

313 ± 16 


J only 

J0857060-055418 

MCG-01-23-010 NEDOl 

6dE 

0.0454 

313 ± 16 


J only 

J1153082-323357 

ES0440-G025 

6dE 

0.0271 

313 ±22 



J1343331-381033 

ES0325-G004 

6dE 

0.0339 

356 ± 16 

SNL-0 

H only 

J1352252-345600 

ES0384-G007 

6dE 

0.0382 

356 ±21 



J1352533-282921 

NGC 5328 

6dE 

0.0158 

335 ± 8 



J1450280+110649 

2MASX J14502809+1106497 

SDSS 

0.0533 

323 ± 6 



J1510556-112847 

NGC 5872 

6dE 

0.0249 

329 ± 14 



J1511314+071507 

CGCG 049-033 

SDSS 

0.0446 

334 ± 6 



J1557496+161836 

NGC 6023 

SDSS 

0.0370 

328 ± 5 



J1601361+122136 

CGCG 079-022 

SDSS 

0.0356 

325 ± 5 



J1613129+174828 

2MASX J16131301+1748283 

SDSS 

0.0372 

316 ± 7 



J1832293-601726 

2MASX J18322937-6017262 

6dE 

0.0517 

358 ± 46 



J1916325M01233 

NGC 6768 

6dE 

0.0186 

310 ±13 



J1958358-523509 

EAIRALL 0877 

6dE 

0.0509 

307 ± 24 



J2100286M22523 

ESO286-G022 

6dE 

0.0312 

356 ± 18 

SNL-1 


J2129374-211144 

IC 1386 

6dE 

0.0356 

319 ±14 



J2318463-102357 

IC 1479 

6dE 

0.0317 

317 ±19 



J2324259+135820 

NGC 7651 

SDSS 

0.0435 

322 ± 5 




and H bands. Hence, we calculated that a survey of only ~30-40 
galaxies is needed to yield a meaningful sample of 3-4 clean IMF- 
constraining low -2 lens systems, and compiled a target sample ac¬ 
cordingly. 

2.2.1 Sample 

To maximize the number of detectable lensed sources, the candi¬ 
date sample is explicitly biased towards the most massive galax¬ 
ies. Targets were selected from the Six-degree Field (6dF) Peculiar 
Velocity Survey (Campbell et al. 2014), which covers the whole 
southern sky with |6| > 10°, and from the Sloan Digital Sky Sur¬ 
vey (SDSS, Data Release 7, Abazajian et al. 2009) in the north, 
with R. A. and declination cuts suitable for observation in the April- 
September period from the VLT. 

The primary selection was on stellar velocity dispersion, 
withagdF >300kms“^ or ctsdss >310kms“^ (the difference ac¬ 
counts for the different fibre sizes in the two surveys, 6.7 arcsec 
vs 3.0 arcsec, given typical velocity dispersion gradients). Camp¬ 
bell et al. measured velocity dispersions only for 6dF galaxies with 
2 <0.055; we apply the same redshift limit on the SDSS candi¬ 
dates. A signal-to-noise limit (15 A~^) was imposed to exclude 
galaxies with large errors on velocity dispersion. These initial se¬ 
lection criteria produced a list of ~130 candidates. 

In clusters or high-mass galaxy groups, the large dark-matter 
component can push the lensing critical curve to large radius. 
Such cases, in which lensing measures the cluster-scale halo, rather 
than the galaxy mass, are much less suitable for constraining the 
stellar IMF. We used the NASA Extragalactic Database, together 
with compilations of group velocity dispersion (o-grp) data (e.g. 
Struble & Rood 1999; Ragone et al. 2006; Crook et al. 2007; 
Coziol et al. 2009) to exclude central galaxies of groups with 


Cgrp > 400kms“^). Additionally, we removed non-central galax¬ 
ies that lie within (or projected close to) very massive clusters, 
to improve robustness of the lensing masses. Targets from 6dF at 
low galactic latitude (|6| < 20°) were also excluded, since clus¬ 
ter/group catalogues are less complete in these regions. 

All available survey images of the targets were inspected to 
remove galaxies with double nuclei or very close companions, 
which might have influenced the measured velocity dispersion (e.g. 
Bernardi et al. 2006). The SDSS and 6dF survey spectra were in¬ 
spected to remove galaxies with strong emission lines, since the 
presence of young stars, or continuum from an active galactic nu¬ 
cleus, would hinder the eventual inference on the IMF. 

Of the 38 galaxies reaching this stage in the selection process, 
only four have been imaged with HST. One of these is the known 
lens ES0325-G004, which was retained in the sample (for H-band 
observation only) as a method verification target. Another object 
with HST data is NGC 5532, where a faint arc-like feature has been 
reported from shallow ultra-violet imaging (Hardcastle et al. 2005), 
but has not been confirmed with spectroscopy; we retained this can¬ 
didate in the target list (although it was not ultimately observed). 
For the final two galaxies (NGC 6768 and ES0508-IG045), the 
HST images do not suggest the presence of any lensed background 
galaxies; these targets were removed from the SINFONI sample. 
The final target list hence comprises 36 lens candidates. 

2.2.2 Observations and data reduction 

The SINFONI observations were made in service mode during 
ESO Period 93 (April-September 2014). A complete observa¬ 
tion for each target consists of two pointings in J band and two 
in H band, with a single 600-sec exposure in each pointing, 
per band. Hence the total exposure time for each galaxy is 40 
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pixel pixel 

Figure 1. Sections from the “detection images” showing detected background line emission in the four galaxies described in Section|^ The detection image is 
a wavelength-calibrated, noise-normalised unwrapped cube, i.e. with the IFU slices placed alongside each other (the slice edges are faintly visible in SNL-2). 
The two pointings have been subtracted, resulting in positive and negative signals, and processed to suppress signal from the lens galaxy and residual sky, 
highlighting the background emission line features. 


minutes (rather than the 60 minutes assumed in the the survey 
yield estimate). Using the largest pixel scale available with SIN- 
FONI, each pointing covers a field-of-view 8x8arcsec^, sampled 
with 0.250x0.125 arcsec^ pixels. The two pointings were offset 
hy 2.3 arcsec, always aligned NE-SW. The spectral resolution is 
~2000 in J and ~3000 in H. At the close of Period 93, a total of 
27 candidate galaxies had heen observed, from the original sample 
of 36. Two candidates were observed only in J band due to techni¬ 
cal problems, while the known lens ES0325-G004 was targetted 
for observation only in H. The complete list of observed targets is 
summarized in Table [T] 

Each pair of offset exposures was reduced using standard 
methods from the ESO SINFONI pipeline, to produce an “un¬ 
wrapped cube”, in which all spatial slices are rebinned onto a cali¬ 
brated wavelength scale, but formatted as a two-dimensional spec¬ 
trum, containing positive and negative traces of the target, together 
with residual airglow signal. No correction for telluric absorption 
was made. Further processing steps specific to SNELLS were ap¬ 
plied to remove the residual sky and the strong continuum from the 
elliptical galaxy, leaving any faint background emission lines un¬ 
touched. The galaxy signal was removed by modelling the source as 
having a spatially-invariable spectrum (derived from the bright cen¬ 
tral region), modulated by a low-order continuum correction fitted 
to each spatial pixel. The median residual sky was constructed and 
subtracted for each spatial slice. After masking bad pixels (based 
on a Laplace transform filter) and the extreme ends of each slice, 
all pixels were normalised by the marginal rms function computed 


in both the spatial and the spectral direction, to yield an approx¬ 
imation to a signal-to-noise image, which was then smoothed to 
highlight features on the scale expected for emission-line features. 
We refer to this product, still in unwrapped cube format, as the “de¬ 
tection image”. For most galaxies this is essentially a blank frame 
with fairly uniform noise. 

Candidate lenses were identified through visual inspection of 
the detection images to locate background emission linefl The 
characteristic features are peaks on the appropriate spatial and spec¬ 
tral scale, usually present both in positive and negative (i.e. in both 
two spatial pointings), extending across multiple image slices, and 
not coincident with the galaxy nucleus. Clear signals of this form 
were found in four cases, as shown in Figure [T] We present details 
of these systems in the next section. 


3 IDENTIFIED LENSES 

After identifying emission lines in the detection image, we extract 
narrow-band images centred on the brightest line, to establish the 
spatial form of the system, and compute spectra for the bright¬ 
est images to determine the source redshifts. Figure shows the 

^ Although a well-designed automated detection algorithm might be able 
to locate fainter sources, such systems would probably require dedicated 
deep follow-up observations to confirm them as lenses, and hence lack one 
of the advantages of the SNELLS approach. 
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Figure 2. Source and lens images, and extracted spectra, for the SNELLS lens candidates. For each source, the first panel shows a narrow-band, continuum- 
subtracted image extracted around the brightest line (Ha for SNL-1, [O III] 5007 A in all others). The lens galaxy centre is mai'ked by the red cross; the compass 
shows the field orientation. The second panel shows the broad-band image as greyscale, with contours indicating the narrow-band flux, and source images 
identified. The third panel shows spectra extracted from the source images. These spectra have been normalised by the errors to suppress sky subtraction 
residuals, hence they are strictly closer to signal-to-noise spectra. The strongest emission lines are indicated. For SNL-O and SNL-3, the data used are the 
original SNELLS observations (two offset exposures). For SNL-2, two repeated observations (four exposures) were combined. For SNL-1, we show the 
deeper follow-up data obtained in director’s discretionary time. 
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Figure 3. Images of the three confirmed SNELLS lenses SNL-0,1,2, and the candidate SNL-3. The upper row shows optical imaging for each galaxy, from 
SDSS where available, and otherwise from SuperCOSMOS digitized Schmidt plates (Hambly et al. 2001). The red outline indicates the '^lOx 10arcsec^ 
footprint of the SNELLS obseiwations. The lower row shows the 2MASS J-band images (Jarrett et al. 2000). 



Figure 4. Optical spectra of the three SNELLS lenses SNL-0,1,2, and the candidate SNL-3, from the 6dF Galaxy Survey (Jones et al. 2004, 2009). The signal- 
to-noise ratio in each spectrum is 20-30 A~^. For SNL-0 we also show the high-quality VIMOS spectrum analysed in SL13, with signal-to-noise ~400 
(red). Grey sections indicate spectral regions affected by sky lines or atmospheric absorption. Green dotted lines indicate the potential location of Ha, H/3, 
[O III] and [N II] emission lines. Weak [N II] emission is seen in SNL-3, but there is no evidence for any nebular emission in SNL-0,1,2. 
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SNELLS data for each galaxy, in the form of (a) the continuum- 
subtracted narrow-band image centred on the brightest detected line 
from the source, (b) a broadband image showing the lens, with the 
narrow-band contours overlain, and (c) the extracted spectra from 
the brightest images of the source. 

In three of the four candidates (one of which is ES0325- 
G004), we recover multiple emission-line sources with identical 
wavelengths; the interpretation of these systems as strong gravita¬ 
tional lenses is compelling. In the fourth case, only a single source 
image is observed. 

For additional context, we show the available optical/IR imag¬ 
ing for the lens candidates in Figure[^and the optical spectra from 
6 dF in Figure |4] The following sections describe each candidate 
lens in turn. 

3.1 J1343331-381033 (ES0325-G004, SNL-0) 

This system is the previously-known low-redshift lensing ellipti¬ 
cal (Smith et al. 2005; SL13), which we assign the name SNL-0 
for the purposes of this paper. Satisfying our SNELLS selection 
criteria, this galaxy was observed only in the H band, in order to 
test methods for data processing and source detection, in a case 
where the presence of lensed emission lines is known in advance. 
SNL-0 is the central galaxy of a poor group (Abell S0740), and 
has boxy outer isophotes but no disturbances in its inner regions. 
Compared to the distribution of sizes and luminosities of high-a 
galaxies, SNL-0 has a large effective radius ( 6.6 kpc, compared to 
median 5 kpc for CT>300kms“^ galaxies from Campbell et al. 
2014) and large total luminosity (4x compared to median 

2 xlO“L© .). The 6 dF optical spectrum (supported by a high-S/N 
VIMOS spectrum from SL13) shows an apparently passive stel¬ 
lar population, with weak Balmer absorption lines and no traces of 
emission. 

The SNELLS data for SNL-0 are shown in the top row of Fig- 
ure|2]. The narrow-band emission-line image clearly recovers most 
of the structure of the arc system as seen in HST images, though 
some fidelity is lost from combining only two offset frames. Our 
observations recover the [Olll] doublet and weak H/I line, as previ¬ 
ously seen in long-slit data from SL13. No other lines are detected. 

3.2 J2100286-422523 (ESO286-G022, SNL-1) 

This is the first newly-discovered lens system from SNELLS (Fig- 
ure|2 row 2). With an effective radius of only 2 kpc, SNL-1 is three 
times smaller than SNL-0, and among the most compact objects 
among the ct > 300kms“^ galaxies from Campbell et al. (2014). 
The total luminosity of SNL-1 is accordingly smaller (~10^^ Lq)- 
Due to its southerly declination, SNL-1 has not been covered by 
any modern public optical imaging survey; the SuperCOSMOS and 
2MASS images show an elongated but otherwise featureless mor¬ 
phology (Figurej^. The 6 dF spectrum suggests an old passive stel¬ 
lar population, with weak H 7 and H/1 absorption, and no emission 
detectable at Ha or elsewhere (Figure|4ll. There are no other galax¬ 
ies or groups close to SNL-1; in fact, this object is included in the 
catalogue of unusually isolated galaxies compiled by Karachent¬ 
seva et al. ( 2010 ). 

The original SNELLS data revealed a background emission 
line in the J-band, forming a source at 3.7 arcsec from the target 
centre and a suspected faint counter image on the opposite side of 
the target, at 1.1 arcsec. Deeper J-band follow-up observations of 
SNL-1 were obtained in a director’s discretionary time (DDT) ob¬ 
servation with SINFONI on 2014 August 12. We observed for a 


total of 60 minutes on-source, with equal time on fully-offset back¬ 
ground areas to provide a cleaner sky subtraction. Three different 
detector orientations were used, to help in rejecting cosmetic de¬ 
fects from the images. The deep and clean follow-up data confirm 
the strongest detected line as Ha at a redshift of z = 0.926. The 
[N II] 6548,6584 A and [S II] 6717,6731 A lines are also detected 
in the brighter image. The fainter counter-image is unambiguously 
detected in the DDT observation, with multiple lines confirming 
the lensing interpretation beyond reasonable doubt. 

With aiens = 0.031, SNL-1 is now the closest confirmed 
galaxy-scale strong lens. 

3.3 J0141423-073528 (2MASX J01414232-0735281, SNL-2) 

This is our second newly-discovered lens system (Figure|^ row 3). 
SNL-2, with 2 iens = 0.052, is part of a galaxy pair (Figure 13; the 
companion galaxy is ~1 mag fainter, located 7 arcsec to the North. 
SDSS imaging shows that the pair is embedded in a common, 
asymmetric stellar halo, indicating the galaxies are at a common 
distance and suggesting a recent or ongoing tidal interaction. SNL- 
2 itself appears very regular; its effective radius ( 6.0 kpc) is com¬ 
parable to SNL-0, and its total luminosity is only slightly smaller 
(SxlO^^L©). The 6 dF spectrum is consistent with a pure passive 
stellar continuum (Figure |4ll, but the Ha/[Nll] region is contam¬ 
inated by telluric absorption so we cannot rule out weak nebular 
emission as securely as for SNL-0 and SNL-l. 

By chance, two separate SNELLS observations were obtained 
for SNL-2, so we combine the data to improve total deptlfl The 
H-band observation reveals emission in two lines consistent with 
the [O III] 4959,5007 A doublet, at a redshift of 2=1.97. (The 
observer’s-frame wavelength is ~ 1.48 pm, placing the lines at the 
short-wavelength edge of the H band, in a region of poor atmo¬ 
spheric transmission.) The lines are noticeably broad, with FWHM 
~550kms“^. The emission is concentrated in a bright source at 
3.2 arcsec from the galaxy centre, with a fainter second image at 
radius 1.2 arcsec diametrically opposite. The 4959-A line is not 
clearly seen in the fainter image, but we consider the concordant 
5007 -A line, and the symmetric configuration, sufficient to confirm 
this as a lensed counter-image. No counterpart to the emission-line 
source is visible in the SDSS imaging. 

3.4 J0258177-044906 (MCG-01-08-023, SNL-3) 

Finally, we report a further candidate lens system, which cannot 
be modelled sufficiently well to include in our mass estimates 
(Figure 12 row 4). In SNL-3, the lens is an edge-on SO galaxy at 
2 = 0.031, with a prominent dust lane seen in SDSS imaging (Fig¬ 
ure [2l. The 6 dF spectrum shows weak emission at [Nil]6584 A 
from the lens candidate itself (Figure|4j. 

The SNELLS data reveal emission lines consistent with the 
[O III] 4959, 5007 A doublet at a redshift of 2 = 1.95 (as for SNL-2, 
this places the lines in a region of low atmospheric transmission at 
the edge of the H-band). The source is 3.9 arcsec from the candi¬ 
date lens centre and appears to be tangentially extended; however, 
no counter-image can be identified with the present data. No coun¬ 
terpart to the emission-line source is visible in the SDSS imaging. 

There are three possible interpretations for SNL-3: (a) the 
background galaxy is not located in the multiply-imaged region 

® No observatory flux standard was observed for the first observation, so 
the target was returned to the queue. 


© 2015 RAS, MNRAS OOP.m-fTT] 


SNELLS low-redshift lenses and the IML 9 


in the source plane, in which case only a weak constraint on the 
lensing mass can be derived; (b) the source is doubly imaged, with 
a faint inner counter-image that is too faint to detect without deep 
follow-up; (c) an outer counter-image lies beyond the SINFONI ob¬ 
servational footprint, implying a very large lensing mass (we dis¬ 
cuss field-of-view selection bias more generally in Section l4^ . 

Although improved data might eventually distinguish between 
these possibilities, the complexity of the lens galaxy (dust lane, 
emission lines) makes this a poor system for inferring the IMF 
mass factor in any case, and no further observations are planned 
for SNL-3. 


4 FROM LENSING TO IMF CONSTRAINTS 

The observed lensed can be used to constrain the IMF mass excess 
factor, a, defined by 

^ _ T _ / Msin — Mbm \ 1 

Tref V I/Ein / Tref 

where MEin is the total lensing mass within the Einstein radius, 
Mdm is the contribution from dark matter, Leiu is the Einstein 
aperture luminosity, and T^ef is a reference mass-to-light ratio for 
the stellar population assuming a fiducial IMF. 

In this section we discuss the measurements and assumptions 
made for each quantity above, and combine them to yield con¬ 
straints on a. The lensing mass and the associated luminosity are 
strongly constrained by the data (Section EU, while Mom and 
Tref involve assumptions which are either not determined at all 
from the observations (dark matter) or only partly informed by data 
(stellar population model). Our default treatment of these quanti¬ 
ties is described in Section |4^ In Section 143] we test alternative 
assumptions and discuss other possible systematic effects on our 
results. Finally, in Section I4A1 we discuss a selection bias due to 
the limited field of view of the IFU observations. 


4.1 Lensing M/L 


In this section we report the derivation of the total projected mass 
and projected luminosity within the Einstein radius of each lens. 

Eor SNL-0, we simply adopt the angular Einstein radius, 
f^Ein = 2.85 arcsec, as determined by Smith et al. (2005) from fit¬ 
ting HST imaging using a singular isothermal sphere mass model, 
with a linear external shear term. 

Eor SNL-1 and SNL-2, we approximate i?Ein by half the im¬ 
age separation, and compute the associated lensing mass using the 
symmetric lens-mass formula, 


MEin = = 1.23 X 10' 

4G Dla 


'Mo {^j 


D 


100 Mpc 


(e.g. Refsdal & Surdej 1994). The effective distance D is calcu¬ 
lated from the angular diameter distances D\ (observer to lens), Dg 
(observer to source) and D\s (lens to source). It is implicit in this 
approach that the lensing potential is assumed to be fixed on the 
centre of light; this is clearly justified if the lensing mass is indeed 
dominated by stars. 

Using the gravlens code (Keeton 2001), we have verified 
that the mass projected within a circular aperture, with radius given 
by the half-separation of the image pair, is well reproduced by ap¬ 
plying the symmetric approximation, even when the image con¬ 
figuration is quite far from symmetric, as in SNL-1 and SNL-2. 





a 


Figure 5. IMF mass factor constraints for the three SNELLS lens galax¬ 
ies. In the upper panel, the heavy solid curves show the probability disti'i- 
bution for a in each galaxy, adopting our default assumptions (old stellar 
population, EAGLE dark matter contributions). The second panel shows 
the results under the extreme assumption of negligible dark matter within 
the Einstein radius. The third panel shows the results obtained when fit¬ 
ting Tref from the 6dF spectra, rather than assuming the populations are 
old. The heavy/light bar below shows the \cr/2a intervals for the mean of 
the three galaxies. For comparison, the yellow shading indicates the ±lo¬ 
in a expected at the mean velocity dispersion of our sample, 345 km s“ ^, 
based on the trend reported by Cappellaii et al. (2013) (slightly extrapolated 
since their sample contains no very high-o- galaxies). The a values for sev¬ 
eral candidate IMFs are shown with grey bars. Here, “heavyweight” refers 
to q; = 2, corresponding to the heaviest results from Treu et al. (2010) and 
Conroy & van Dokkum (2012b). 
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Table 2. Observed quantities, con'ections and derived parameters for the three confirmed SNELLS lenses. 


Name 

SNL-0 

SNL-1 

SNL-2 

notes 

Target 

J1343331-381033 

12100286^22523 

10141423-073528 


NED ID 

ES0325-G004 

ESO286-G022 

2MASX101414232-0735281 


■^lens 

0.034 

0.031 

0.052 

lens redshift 

^^src 

2.141 

0.926 

1.969 

source redshift 

Dl [Mpc] 

149 

135 

230 

luminosity distance to lens 

Da [Mpc] 

139 

127 

208 

angular diameter distance to lens 

/ = 

1.027 

1.043 

1.044 

lensing geometry factor 

D [Mpc] = fDA 

143 

132 

217 

effective distance to lens 

Scrit [1O®M0PC“2] 

12.31 

13.05 

7.35 

lensing critical surface density 

o-6dF [kms“l] 

356±16 

356±18 

320±18 

aperture corrected to i^eff /8 

[Mg/Ee] 

0.33±0.05 

0.31 ±0.05 

0.38L0.06 

full-spectrum fit to 6dF 

i?Eft [arcsec] (J) 

9.8 

3.3 

5.9 

PSF-corrected 

fEff [kpc] (J) 

6.6 

2.0 

6.0 

PSF-corrected 

Txot [lO^°L0,j] 

40.9 

9.8 

25.5 


J?Ein [arcsec] 

2.85 

2.38 

2.21 

0.5 X image-separation for SNL-1,2 

rEin [kpc] 

1.92 

1.47 

2.23 


fcorr 

1.00 

0.96 

0.88 

correction factor for shear / asymmetry 

MEin [1O1OM0] 

14.3±0.7 

8.9L0.4 

11.5LI.2 


'^Ein 

12.36 

12.80 

13.53 

PSF-corrected 

Aj 

0.04 

0.03 

0.02 

Schlafly & Finkbeiner (2011) 

kj 

-0.02 

-0.02 

-0.03 

bandshifting correction (via EzGal) 

LEin [IQIOLq.j] 

7.8±0.2 

4.3L0.1 

6.2±0.2 

assumes Mq^j = 3.71 as in Maraston (2005) 

M/L [solar units, J-band] 

1.82±0.10 

2.08L0.12 

1.87L0.22 

total, i.e. includes DM 


For isothermal and de Vaucouleurs profile mass models with el- 
lipticities up to e = 0.25, we find a scatter of < 3 per cent, and a 
bias < 1.5 per cent. External shear introduces a comparable level 
of noise. Including all sources of error, we estimate uncertainties of 
~5 per cent on the lensing masses. 

For both SNL-1 and SNL-2, we apply small corrections to the 
symmetric lens approximation. In SNL-I, the light distribution has 
clear ellipticity and the images are aligned along the minor axis. 
In this configuration, if the lensing mass is distributed with similar 
shape and orientation as the light (as must be the case if stellar mass 
dominates within i?Bin), then the lensing mass is overestimated 
by applying the symmetric lens approximation. For isothermal and 
de Vaucouleurs profiles with ellipticities e = 0.2-0.3, we find that 
MEin should be revised downwards by 4 per cent for SNL-1. 

For SNL-2, there is an additional complication due to the close 
neighbouring galaxy, which must induce some distortion to the ob¬ 
served lens configuration. We explore the probable influence of the 
companion using gravlens. We build a two-component model, 
with singular isothermal sphere mass profiles fixed at the coordi¬ 
nates of SNL-2 and the neighbouring galaxy 7 arcsec north. With 
increasing mass contributed by the neighbour, the source-plane (un- 
lensed) location of the emission line object must move northwards, 
and the normalization of the primary lens must be reduced, in or¬ 
der to maintain the same observed image configuration. The true 
mass projected within the “naive” i?Ein (the half-separation of the 


images) is also smaller if we allow for substantial mass in the com¬ 
panion. Hence the effect of neglecting the external shear would be 
to overestimate the lensing mass of SNL-2. If the companion has 
10-50 per cent of the mass of the primary lens (as suggested by 
the flux ratio), then the lensing mass for SNL-2 should be revised 
downwards by 8-16 per cent. We adopt a correction of 12 per cent 
and an increased lensing mass error of 10 per cent to account for 
greater uncertainty in the external shear for this system. 

With these corrections applied, the Einstein masses for SNL-0, 
SNL-1 and SNL-2 are 14.3, 8.9 and 11.5 xlO^'^ Mq, respectively, 
with adopted errors of 5, 5, and 10 per cent. 

To derive the Einstein aperture luminosity Leiii on a common 
basis for all three lenses (two of which - SNL-0 and SNL-1 - are 
beyond the coverage of SDSS), we use J-band photometry from 
2MASS (larrett et al. 2000). The J band is preferred over H or K 
for this purpose, because the fluxes will ultimately be compared 
to stellar population models, which show better agreement among 
themselves at 1 than at longer wavelengths (see e.g. fig 4 of Man- 
cone & Gonzalez 2012). 

Aperture fluxes were measured directly from the image data 
within the Einstein radius of each lens. A correction for the point- 
spread function (PSF) was derived by fitting a convolved Sersic 
model using galfit (Peng et al. 2010). The convolution ker¬ 
nel was derived from stars in the same image frame as the target 
galaxy. The PSF corrections are ^0.25 mag, and we estimate them 
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to be accurate to 0.02 mag. This uncertainty is added in quadrature 
with photometric errors tabulated in the 2 MASS database which, 
for such bright galaxies, are dominated by systematic calibration 
errors of ~0.02-0.03 mag. Small corrections are applied for band- 
shifting (“k-correction”), based on estimates from various mod¬ 
els implemented in EzGal (Mancone & Gonzalez 2012), and for 
galactic extinction, from Schlafly & Finkbeiner (2011). To convert 
from flux to luminosity, we use the Hubble-flow luminosity dis¬ 
tance in the adopted cosmology, and the absolute magnitude of the 
Sun Mj = 3.71, as assumed in the Maraston (2005) models that we 
use to estimate Tref below. The resulting Lbiii are 7.9, 4.3 and 

6.1 X lO^°Z /0 for SNL-0, SNL-1, and SNL-2 respectively, with to¬ 
tal estimated errors of 3^ per cent. The total J-band mass-to-light 
ratios (including dark matter) for the three galaxies are therefore 
1.8±0.1, 2.1±0.1 and 1.9±0.2. 


4.2 IMF normalisation factor 

In this section we describe the derivation of the IMF factor, a, using 
our “default” assumptions both for the dark matter content and for 
the stellar population model. 

The dark-matter corrections within iZEin are expected to be 
small, since the lensing measurement probes the baryon-dominated 
inner part of the galaxy. In SL13, we estimated the correction for 
SNL-0 based on halo properties in pure dark-matter simulations. 
We assumed Navarro, Frenk & White (1996) halo profiles, with 
the mass-concentration relation of Neto et al. (2007), and impos¬ 
ing limits on the halo mass according to the observed velocity dis¬ 
persion of the surrounding group. The limitations of the approach 
adopted in SL13 are (i) that the simulations do not trace any re¬ 
sponse of the halo the to the growing presence of a dense stellar 
system at its centre and (ii) that the calculation assumes the scatter 
in the mass-concentration relation is uncorrelated with the likeli¬ 
hood of hosting a giant elliptical galax} 0 . 

For the present work, we instead use dark matter masses es¬ 
timated directly from the EAGLifl cosmological hydrodynamical 
simulation (Schaye et al. 2015). This simulation includes state- 
of-the-art baryonic physics prescriptions and is able to reproduce 
a wide range of observables for the local galaxy population. En¬ 
closed (i.e. 3D) density profiles from EAGLE have been presented 
as a function of halo mass by Schaller et al. (2015), who showed 
that the inclusion of baryons significantly affects the central re¬ 
gions. To ensure that the corrections are as well matched as possi¬ 
ble, Schaller kindly recomputed projected profiles for galaxies for 
galaxies that match the selection of lens candidates in SNELLS. 
Specifically, we selected simulated halo-central galaxies with stel¬ 
lar mass IO^'^-IO^^Mq and projected stellar velocity dispersions 
a > 275 but excluding halos with total mass greater than 

ltL*MQ (to reject massive groups and clusters). With this selec¬ 
tion, we find that the projected dark matter masses from EAGLE, 
interpolated to match the Einstein radii for our lenses, are remark¬ 
ably similar among galaxies. The mean EAGLE Mdm measured 
at the Einstein radius for SNL-0, SNL-1 and SNL-2 is 2.3, 1.4 
and 3.Ox 10^°M© respectively, with a halo-to-halo scatter of only 


® Neto et al. (2007) show that much of the scatter around the 
concentration-mass relation is attributable to variation in halo formation 
redshift. If ellipticals are preferentially found in the earliest-forming halos, 
the SL13 method would underestimate the dark-matter content on average, 
and overestimate the galaxy-to-galaxy scatter. 

^ Evolution and Assembly of GaLaxies and their Environments. 


17 per cent in each case. Expressed as a fraction of the total lensing 
mass, the dark-matter correction is I 6 I 2 SNL-0 and 

SNL-1, and 26^4 per cent for SNL-2. 

Other possible non-stellar mass contributions within the Ein¬ 
stein radius, e.g. gas and central black holes, are neglected. We 
implicitly assume that all gas lost during stellar evolution has ei¬ 
ther been heated and dispersed into a low-density halo or else re¬ 
cycled into subsequent generations of stars, and hence does not re¬ 
main as an additional mass component to be subtracted. The aver¬ 
age central black-hole mass for a = 300 kms“^ early-type galaxies 
is ~L4x 10® M© (McConnell et al. 2011), an order of magnitude 
smaller than the dark-matter contribution. 

The reference stellar mass-to-light ratio, Tref, is the ratio ex¬ 
pected for each galaxy assuming a universal fiducial IME, here 
taken as that of Kroupa (2001). The T^ef depends mainly on galaxy 
“age”, in a luminosity-weighted sense, with a weak effect from 
metallicity. In SL 13, we fitted to Lick indices measured on a high- 
S/N spectrum of SNL-0, to determine age and metallicity con¬ 
straints, from which we derived limits on Tref. For the enlarged 
sample of three lenses, we do not have equivalently high-quality 
spectra for all of the galaxies. For the “default” results of this pa¬ 
per, we instead proceed by assuming all three galaxies are “old”, 
which is supported by analysis of larger galaxy samples (e.g. Mc- 
Dermid et al. 2015). Specifically, we adopt Tref values appropriate 
for a mean age of lOGyr, and allow a standard deviation of 1 Gyr. 
The corresponding 2 ct range for formation redshift is 1 < 2 < 4. At 
solar metallicity, the predicted distribution for (rest-frame J-band) 
Tref from Maraston (2005) models has mean 1.47 and standard de¬ 
viation 0.12. A factor of two in metallicity yields only a 3 per cent 
change in Tref in J. 

Finally, we combine the above quantities and their probabil¬ 
ity distributions to arrive at estimates for a, relative to the Kroupa 
IMF: a = L04±0.11, L20±0.13 and 0.94±0.17, for SNL-0, SNL- 
1 and SNL-2, respectively. The larger error for SNL-2 reflects the 
larger dark-matter correction necessary due to its larger distance, 
and also the more uncertain shear due to the companion galaxy. 
The result for SNL-0 is consistent with that obtained in SL13 
(L04±0.15), despite several changes to the assumptions, input data 
and methodology. 

The principal result of this paper, then, is that all three known 
low-redshift early-type lens galaxies have mass-to-light ratios con¬ 
sistent with a Kroupa IMF. Each individually is incompatible 
with a Salpeter (a= 1.55) or heavier IME, at the >99.5 per cent 
level, with respect to the random uncertainties (we discuss sys¬ 
tematic errors below). The average a for these three very massive 
(ct > 300kms“^) early-type galaxies is (a) = L06±0.08. This re¬ 
sult is summarized in Part (a) of Table and in the upper panel of 
Figure|5] 


4.3 Robustness and systematics 

According to the default assumptions made in the previous section, 
our results favour lightweight IMFs in all three confirmed SNELLS 
lenses. In this section we explore the sensitivity of the result to 
modifications in our methodology and inputs. 

The lensing mass and corresponding luminosity, both pro¬ 
jected quantities within the nominal Einstein radius, are the most 
robustly determined properties of the systems. Neither quantity 
could be systematically in error at the ~50per cent level neces¬ 
sary to derive a Salpeter IME. Instead, the main uncertainties arise 
in correcting the total lensing mass to a stellar mass, and in ref- 
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Table 3. Conversion from total leasing mass-to-light ratio to IMF mass factor a. We summarize the limits and/or assumptions made for the projected dark- 
matter mass (within the Einstein radius), Mqm. and for the reference mass-to-light value, Tj-ef, and report the resulting constraints on a. Part (a) shows the 
assumptions and results for our “default” treatment (see Section l4~^ . while subsequent parts show how the results are affected by changes to our assumptions 
and methods (Section|43)- The final column reports the mean a for the three lenses; including a systematic error component of 10 per cent. 



Assumptions 

Quantity 

SNL-0 

SNL-1 

SNL-2 

{«> 

(a) 

Eagle DM 

Mdm 

Q QQ”t"0.38 
z.zo_o 32 

1 3Q+0-24 

l.o»_q 20 

^•^®-0.42 



Assumed old (M05 models) 

Tref 

I.47±0.12 

1.47±0.12 

1.47±0.12 



(default result) 

a = T/Tref 

1.04±0.11 

1.20±0.13 

0.94±0.17 

1.06±0.08±0.10 

(b) 

SL13 dark matter method 

Mdm 


1 1 ^+0-85 
^•^‘^-0.49 

9 C7+1-90 

-1.09 




« = T/Tref 


1 /^q+0.15 
^•^■^-0.19 

0 99+0-20 
^■^^-0.27 

1.10+°;“±0.10 

(c) 

No dark matter 

« = T/Tref 

^•^^-0.11 


1 97+0.19 
^•^'-0.17 

1.31±0.09±0.10 

(d) 

Eit ages from 6dE (CvD12 models) 

Tref 

1.27 ±0.23 

1.46 ± 0.18 

1.40 ±0.27 




a = T/Tref 

1 9C+0.29 
^•^•^-0.21 

1 ^c+019 

16 

i.uj o 23 

i.i9/;°;i®±o.io 

(e) 

Allow for FoV selection bias 





1.10±0.08±0.10 


erencing the stellar mass-to-light ratio to the value expected for a 
fiducial IMF. 

The dark-matter mass contribution is an uncertain input to 
the calculation, as it is not estimated directly from the data. How¬ 
ever, since the correction factor is small, plausible variations in the 
dark matter treatment make only small impact on a. To illustrate 
this, we have repeated the calculation using estimates based on 
the concentration-mass relation in pure dark-matter simulations, 
as in SL13. In this case, the mean projected dark-matter masses are 
~20per cent smaller than in the default model The SL13 scheme 
also yields much larger scatter in dark-matter content, and hence 
a slightly increased error in a. Adopting this approach in place of 
the EAGLE scheme leads to (a) = l.lOlg jj (Tablej^ part b). Even 
the making the limiting assumption of negligible dark matter within 
the Einstein radius, the average derived a would increase only to 
(q) = 1.31±0.09 (see Tablej^ part c, and the second panel of Fig¬ 
ure [Sjl. 

A second input that is not yet derived from the data is the 
assumption that all of the lens galaxies have uniformly old stellar 
populations. Younger galaxies have lower mass-to-light ratios for a 
given IMF, and hence assuming more recent star formation would 
yield a larger estimate of a. In order to favour a Salpeter IMF, we 
would need T^ef ~ LO, implying ages ~6Gyr, on average, for the 
three galaxies. While late-epoch star-formation is known to occur 
in some early-type galaxies (evidenced by emission lines, strong 
Balmer absorption, central dust lanes or ring/spiral structures) this 
is rarer among those with very high velocity dispersion. There is 
no detectable emission in the 6dF spectra for SNL-0, SNL-1 or 
SNL-2, nor do any of them visually show strong Balmer absorp¬ 
tion lines suggestive of recent star formation, i.e. bursts within the 
past ~1 Gyr. Intermediate stellar ages, however, cannot be so easily 
excluded. 

As an alternative to assuming old ages, we have also derived 
probability distributions for T^ef based on fitting the 6dF spec¬ 
tra over the 4000-5400 A wavelength range, using the latest ver¬ 
sion of the Conroy & van Dokkum (2012a, CvD12a) stellar popu¬ 
lation models. Given the limited signal-to-noise (20-30 A“^), we 
fit a single burst stellar population with a restricted set of element 
abundances. Confidence intervals on Tret were derived through a 
Markov Chain Monte Carlo sampling method. Note that the 6dF 
spectra were obtained through a fibre of diameter 6.7 arcsec, and 
hence include some light from beyond the Einstein radius, but 


the mass-to-light ratio from the integrated aperture is unlikely to 
be significantly biased by this mismatch. From the spectral fit¬ 
ting approach fits we derive J-band Tref, assuming a Kroupa IMF, 
of 1.27±0.23, 1.46±0.18 and 1.40±0.27, for SNL-0, SNL-1 and 
SNL-2, respectively. Note that the errors are larger by a factor of 
two than the assumed range of T^ef in our default treatment, im¬ 
plying that the 6dF data quality is not sufficient to confirm that the 
galaxies are ‘old’, although they are certainly consistent with being 
so (the fitted ages are 10±3, 13±3 and 14±4 Gyr). 

The possibility of younger ages leads to a tail in the probability 
distribution towards larger a (third panel of Figure^. The resulting 
mean a over the three galaxies is (a) = l.lOjjg j® (Tabled part 
d). (We account consistently for the different J-band solar absolute 
magnitude used in the CvD12a models, Mj = 3.67, as compared to 
the M05 models used above.) For all galaxies, the spectrum-fitting 
Tref is consistent with the previous assumptions, and hence yields 
a values similar to the default result, though with larger uncertainty. 
We note that the assumption of a single-burst star-formation history 
would lead to a slight overestimate of a using this method, if a 
galaxy actually harbours a more complex population: the V-band 
derived single-burst-equivalent age is more sensitive than the mass- 
to-light ratio to young stellar component, hence at fixed “age”, a 
composite population has larger Tref than a single-burst. 

Our results are subject to various sources of systematic errors 
in stellar population models. Mancone & Gonzalez (2012) report 
a total spread of -^15 per cent between different stellar population 
models in the J-band mass-to-light ratio; of course the various mod¬ 
els share several common ingredients and fundamental calibrations, 
so this does not necessarily reflect the full range of systematic un¬ 
certainty. One specific concern is that most stellar population mod¬ 
els (including Maraston 2005 as used in our default calculation) do 
not fully take into account the effects of non-solar element abun¬ 
dance patterns. Broadband colours can be affected directly by ab¬ 
sorption features in the stellar atmosphere (especially in the blue) 
and by changes in the isochrones due to different opacities in the 
stellar interior. Coelho et al. (2007) presented models taking both 
effects into account consistently, for a range of iron abundance 
Fe/H and a element enhancement. Comparing their 10 Gyr SSP 
models we find that at fixed total metallicity (as assumed in M05), 
a 0.3 dex increase in [a/Fe] yields a -^Tper cent increase in J-band 
Tref- Accounting for this effect would lead to a smaller derived a. 

Finally, we note some more general sources of systematic er- 
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Figure 6. Field-of-view completeness fraction, as a function of both Ein¬ 
stein radius, -Rsin. and IMF mass factor, a. The latter assumes mass profiles 
informed by the 6dF luminosity profiles fore > 300kms“^ galaxies. The 
filled symbols at the bottom of the panels show f?Bin and a for SNL-0,1,2. 
Red dashed lines indicate the 90 and 50 per cent limits. 


ror. Adopting Hq = 67 km Mpc~^ from the Planck cosmology 
(Planck Collaboration 2014), rather than TCq =70kms“^ Mpc“^ 
from WMAP7, would yield smaller a by -^5 per cent. Peculiar ve¬ 
locities of ~500kms“^ would induce shifts of up to 5 per cent 
in a, depending on lens redshift. If the lens galaxies harbour very 
massive central black holes, then a would be a smaller than esti¬ 
mated by our method: As a limiting case, if all lenses are assigned 
Mbh = 10^° Mq, comparable to the most-massive known objects 
(McConnell et al. 2011), (a) is reduced to 0.94. 

Considering the range of values obtained above, we conclude 
that the dominant systematic uncertainty on (a) arises from the 
choice of whether to assume old ages or to fit Tret from 6dF spec¬ 
tra; the two methods yield results that differ by 13 per cent. For 
clarity, we quote the results from both treatments in the abstract. 
The remaining systematics can be absorbed into an additional error 
component of ±0.10, dominated by the variation among different 
stellar population models, but also including the effect of different 
dark-matter assumptions, black hole contributions, abundance ratio 
effects, cosmology etc. No single factor seems likely to account for 
the ~50 per cent discrepancy between SNELLS and studies which 
indicate a heavy IMF in massive galaxies. 


4.4 Field-of-view selection bias 

The systematic errors reviewed in the previous section are uncer¬ 
tainties associated with measurements of a for the detected lenses. 
Here we consider the more subtle effect of a-dependent bias in se¬ 
lection of the lenses. This arises because the SINFONI field of view 
limits the observed area around each target, so that large and 
hence large a, systems may be lost from the sample. In particular, 
for two-image systems like SNL-1 and SNL-2 (more common than 
quadruple images or rings which require closer alignment of lens 
and source), the more distant image may be formed at up to a max¬ 
imum of 2 J^Bin- Although a second image will be formed at small 
radius and hence remain within the field of view, this is the demag- 
nified image, so a large i7Bin system is less likely to be detected 
than a small 7?Bin system. 

To estimate the sample completeness as a function of 7?Bin, 
we used a singular isothermal sphere mass model, with a range 
in normalization, to compute image positions for a uniform grid 
of source locations, relative to the SNELLS observation geometry. 
All systems with i7Bin ^ 2.4arcsec are recovered (i.e. both im¬ 
ages within the field-of-view), after which the completeness drops 
with increasing 7?Bin: the recovered fraction T’(i?Bin) is 90, 50 and 


10per cent at 7?Bin = 2.6, 3.1 and 4.1 arcsec respectively (Figurej^ 
left). 

The connection between i?Bin and a depends on the mass 
profile of the galaxy, and hence estimating the completeness func¬ 
tion for a requires a model for the galaxy-to-galaxy variation in 
the profiles. For this purpose, we use the half-light radius and 
surface brightness information for the 6dF survey galaxies with 
cr > 300kms“^, from Campbell et al. (2014). We assume that the 
inner mass profile is dominated by the stellar mass, following an 
profile with a stellar mass-to-light ratio T = TrefQ:= L5a, 
and allow a constant 15 per cent contribution of dark matter. For a 
given a, the distribution of profiles leads to a distribution of i?Bin, 
from which the completeness function can be derived. As shown in 
the right-hand panel of Figure Plot) fells gradually from unity 
at a = 0.9, to 90, 50 and lOper cent at q= 1.15, 1.70 and 2.70 re¬ 
spectively. 

Combining the completeness F{a) with an assumption for the 
underlying true distribution of a, we can determine the extent of 
the bias in (a). For a case tuned to resemble the dynamical results 
from Cappellari et al., with (a) = 1.60 and scatter 0.3, the mean 
recovered a is 1.52, i.e. a bias of 0.08. For lower input a, where 
the completeness function is shallower, the bias is smaller. For the 
same assumed scatter of 0.3, our recovered (a) = 1.06 is compatible 
with a true underlying mean of 1.10 (Table|3 part e). 

We conclude that although there is a selection bias against 
high-a lenses due to limited field-of-view, this effect is un¬ 
likely to account fully for the small average value obtained from 
SNELLS, compared to results from other studies. Our final esti¬ 
mate for the underlying mean a, correcting for the sample bias, 
is (a) = L10±0.08±0.10 when the lenses are assumed to be old, 
and (a) = 1.231 q ± 0.10 when fitting ages from the 6dE spectra, 

where the second error indicates the remaining systematic compo¬ 
nent. 


5 DISCUSSION 

Figure [7] presents the a estimates for the SNELLS lenses in the 
context of other recent studies, from spectroscopy (CvD12b), dis¬ 
tant lenses (Treu et al. 2010, TIO) and stellar dynamics (Cappellari 
et al. 2013, A3D). The figure is an updated version of that shown 
in SL13. We show a as a function of velocity dispersion and of the 
[Mg/Fe] abundance ratio (a proxy for star-formation timescale), the 
latter being estimated from fits to the 6dF spectra using the CvD12a 
models, as described in Section l43] Both quantities have been con¬ 
sidered in previous work as possible predictors for IMF variations 
(CvD12b; Smith 2012; McDermid et al. 2014). 

Figure |7] confirms that the SNELLS galaxies have lensing- 
derived a factors that are smaller than the average values derived by 
the other studies for galaxies of similar propertie|3- In this section, 
we discuss the comparison with each of the other studies shown in 


® Note that the a-vs-cr relation could in principle be affected by a subtle 
lensing selection bias: lensing is most effective for galaxies with long axis 
aligned towards the observer, maximizing projected mass. While a is not 
biased by this effect (since the projected luminosity is also increased), the 
velocity dispersion for such galaxies would be biased high, relative to a 
random orientation. The [Mg/Fe] abundance ratio is of course not affected 
by such a bias. The high [Mg/Fe] for our galaxies, which lie on the fairly 
tight [Mg/Fe]-(T relation, suggests they do in fact have intrinsically large 
velocity dispersions. 
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the figure, and also comment on compactness as a possible predic¬ 
tor for IMF variation. 


5.1 Comparison with ATLAS3D 

For the comparison with ATLAS3D (A3D) dynamical estimates of 
a, we show the correlation with a as reported by Cappellari et al. 
(2013) and the trend with [Mg/Fe] from McDermid et al. (2014), 
with the estimated intrinsic scatter estimated in those papers. 

The A3D sample does not include any galaxies with 
(j > 300kms“^, but an extrapolation of their trend with a 
would indicate a =1.60 at the mean a of the SNELLS galax¬ 
ies (345kms“^ at i?eff/8, equivalent to (Te = 315kms“^). Tak¬ 
ing into account their estimated intrinsic scatter of 0.084 detfl 
the observed mean of 1.10±0.08 (for the default, “assumed old” 
results, corrected for field-of-view bias, but neglecting system¬ 
atic errors for now) deviates from the A3D relation by A = 2.8(7, 
and Pr(A)= 0.998. Including the systematic error component of 
±10per cent in a reduces the offset of the mean to 2.3 (j; this is a 
lower limit, because some of the systematic error sources are com¬ 
mon to both methods. A test, which does not take account of the 
direction of offsets from the trend, does not detect such a significant 
disagreement. Table |4] summarizes this comparison, together with 
equivalent tests for other a-a relations quoted by Cappellari et al., 
as well as by Posacki et al. (2015) using a combined A3D-I-SLACS 
sample. We conclude that the size of discrepancies from the A3D 
a-(j trend are not especially surprising; however, when taking into 
account that all three galaxies are offset in the same direction, there 
is tension at the 2-3 a level. 

At the average [Mg/Fe] for the SNELLS lenses (0.34), the Mc¬ 
Dermid et al. trend predicts a= 1.50 with intrinsic dispersion of 
0.076 dex. Hence the SNELLS mean a is discrepant from the A3D 
a-[Mg/Fe] relation at a similiar ~2.5 (t level. 

We conclude that although our results significantly favour a 
Milky-Way-like a over a Salpeter-like value for all galaxies, the 
present SNELLS sample is only marginally inconsistent with the 
A3D results, when allowing for the substantial intrinsic scatter 
around their trends. 


5.2 Comparison with SLACS 

At face value, the most surprising discrepancy in Figure [7] is that 
between SNELLS and the results of TIO from SLACS, which is 
also based on strong lensing. In particular, for twelve galaxies with 
(7 > 300 km (after approximate correction to the the same aper¬ 
ture i?eff/8 as used by A3D and CvD12), TIO find (a) « 2.0 with 
no measurable intrinsic scatter. 

Our methods most clearly depart from the approach taken by 
TIO in the treatment of dark matter contributions. TIO fit a spheri¬ 
cal two-component (dark plus stellar) mass model to the measured 
velocity dispersion and Einstein radius. The velocity dispersion de¬ 
pends mainly on the enclosed mass at small radius, while i?Bin con¬ 
strains the projected mass at larger radius. These observables are 
employed together to disentangle the contributions of the two com¬ 
ponents of the mass model. In the TIO analysis, the dark matter is 
assumed to follow the Navarro et al. (1995) density profile, i.e. with 
p(r) r~^ in the region sampled by the data. Posacki et al. (2015) 


® For these comparisons, we assume that the scatter around the A3D trend 
is gaussian in log a. 


show that essentially identical results are obtained using more gen¬ 
eral axisymmetric models constrained by the observed light distri¬ 
bution, velocity dispersion and Einstein-aperture projected mass. 

In SNELLS we have instead chosen to assume a dark matter 
contribution informed by results from the EAGLE hydrodynamical 
simulation. To test whether this difference in method accounts for 
the divergent results obtained, we have used the i?Bin tabulated by 
Auger et al. (2009) to derived the projected dark-matter mass pre¬ 
dicted for the SLACS lenses by our method. For the twelve galaxies 
with the highest velocity dispersions, we find that the dark-matter 
masses predicted from EAGLE are larger, by a median factor of 
2.6 than those recovered by TIO from their two-component mass 
models. Using the EAGLE approach instead would increase the 
median dark-matter fraction (within i?Ein) in the SLACS galax¬ 
ies from 14 per cent to 33 per cent, and reduce the median IMF 
mass factor from 2.0 to 1.6, for these high-(7 galaxies. (Across the 
whole SLACS sample, the median a is reduced from 1.6 to 1.3 
with the imposition of EAGLE dark-matter content.) This reduc¬ 
tion in a is similar to what was obtained by Auger et al. (2010) 
from a re-analysis of the SLACS data, allowing explicitly for halo 
contraction effects. In particular, when modelling contraction fol¬ 
lowing the Gnedin et al. (2004) prescription. Auger et al. recover a 
shallower trend with galaxy mass, reaching a maximum of q « 1.4. 

A second methodological difference between the studies is in 
the method for assigning the reference mass-to-light ratio. Tret: 
TIO use estimates from fits to the spectral energy distribution 
(SED) based on broad-band optical and (sometimes) IR data, de¬ 
scribed by Auger et al. (2009). In fitting the SEDs, however, they 
impose an informative prior on the star-formation histories restrict¬ 
ing the onset of star-formation to early epochs (1 < a < 5, uniform 
in time) and favouring short star-formation timescales (t ^ 1 Gyr). 
In practice these priors force the lens galaxies to have old stellar 
populations, and hence the SLACS Tret is within ~ 10 per cent of 
what would be imposed by our default assumption for SNELLS. 

We have not found a satisfying explanation for the factor-of- 
two disagreement between the lensing-based estimates of a from 
SNELLS and from SLACS for similarly-high-(7 ellipticals. Halo 
contraction could account for part of the differencial but analyses 
of SLACS and other datases appear to favour uncontracted NFW 
halos (e.g. Dutton et al. 2013; Dutton & Treu 2014, but see Grillo 
2012; Sonnenfeld et al. 2012). While it is possible that our dark- 
matter corrections from EAGLE could have been overestimated, 
we note again that even assuming no dark matter at all, the mean a 
from SNELLS is still inconsistent with the average of ^2 derived 
for SLACS at (7 > 300 km s" ^ 


5.3 Comparison with CvD12b 

As emphasised by Smith (2014), there are important differences 
between the A3D and CvD12b results for a, both at the level of 
individual galaxies, and in terms of trends with respect to a and 
[Mg/Fe]. In particular, in CvD12b, a is more strongly correlated 
than with a, while the reverse holds for A3D. For galaxy properties 
similar to those of the SNELLS lenses, CvD12b find larger values 
of a than A3D, with an average of ~2, and hence the apparent dis¬ 
crepancy between CvD12b and SNELLS is (assuming same 
intrinsic scatter as for A3D). Since publication of CvD12b, Conroy 

Since the SLACS analysis incorporates the velocity dispersion, which 
responds to the enclosed mass, it is more sensitive to halo contraction than 
the SNELLS method which uses only projected quantities. 
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Table 4. Compaiison of SNELLS results against the trends reported by Cappellaii et al. (2013) for A3D and Posacki et al. (2015) for a combined A3D and 
SLACS sample. The trends are taken from figure 13 (bottom and top panels) of Cappellari et al. and equations 5 and 6 of Posacki et al. (2015). For each of the 
proposed a-a relations, {log(a)) is the average logo predicted at the mean velocity dispersion of the SNELLS sample and 5 is the intrinsic scatter around 
the trend. The following two columns summarise a test for the three SNELLS observations with respect to the given trend, accounting for measurement 
error in o and a, field-of-view (FoV) bias, and the intrinsic dispersion 5. (We use (5 = 0.081 for intrinsic scatter in the the Posacki et al. trends, based on private 
communication from the authors.) The final two columns are for a comparison of the mean: A is the deviation of measured mean log a from the prediction, 
in units of the error again including field-of-view bias and intrinsic dispersion <5. 


Trend 

(log(a)) 

(5 


Pr(x^) 

A 

Pr(A) 

SNELLS (assumed old; FoV-corrected) 

0.041 ±0.031 






Cappellaii et al. (2013) “prefeiTed” 

0.203±0.048 

0.084 

8.26 

0.9591 

-2.84 

0.9977 

Cappellaii et al. (2013) “all galaxies” 

0.167±0.048 

0.083 

5.18 

0.8406 

-2.22 

0.9869 

Posacki et al. (2015) linear 

0.206±0.047 

0.081 

9.04 

0.9711 

-2.98 

0.9986 

Posacki et al. (2015) quadratic 

0.234±0.047 

0.081 

11.21 

0.9933 

-3.47 

0.9997 




[Mg/Fe] 


Figure 7. The IMF a factor versus u and [Mg/Fe] from SNELLS, compared to previous results from spectroscopy, dynamics and distant lenses. The SNELLS 
results using our default analysis are shown with filled black symbols. Blue open points show the results when deriving from the 6dF spectra; red triangles 
show the extreme assumption of zero dark matter. In the comparison data from Conroy & van Dokkum (2012b), the orange point is M87 and the red point 
is a stack of four Virgo cluster galaxies with high velocity dispersion, from Conroy &. van Dokkum (2010). For ATLAS3D, we show regions corresponding 
to the fitted trends and estimated ilcr intrinsic scatter from Cappellari et al. (2013) and McDermid et al. (2014). The lighter yellow zones coiTespond to 
extrapolations beyond the a and [Mg/Fe] ranges sampled by ATLAS3D. The purple line shows the quadratic fit to a combined ATLAS3D + SLACS sample 
by Posacki et al. (2015). Velocity dispersions have been corrected to an aperture of i?efj/8 throughout. 


& van Dokkum (in preparation) have expanded their observational 
sample, and made numerous improvements to the spectral fitting 
method. These updates have led to somewhat smaller derived a 
factors, on average, for the most massive galaxies. More detailed 
comparisons should await publication of the revised spectroscopic 
constraints. 


5.4 Compactness as a predictor for IMF variations? 

Aside from velocity dispersion and [Mg/Fe], another property of 
galaxies which has been proposed as being related to IMF varia¬ 
tion is “compactness”. Since compact massive galaxies are abun¬ 
dant at high redshifts but rare in the local universe (e.g. Trujillo 
et al. 2006), it has been argued that low-redshift compact galaxies 
are surviving “pristine” relics of the early violent starburst epoch, 
where non-standard IMFs may have been imposed (Lasker et al. 


2013; Chabrier, Hennebelle &. Chariot 2014). In SL13, we argued 
that a dependence on compactness might account for the difference 
between the “lightweight” IMF in SNL-0 and the heavy IMF found 
by Lasker et al. in the 2 = 0.1 compact galaxy “bl9”. 

The addition of the two new SNELLS lenses does not support 
this speculation. Figure[^shows projections of the 6dF Fundamen¬ 
tal Plane from Campbell et al. (2014), and the location of the three 
SNELLS lenses for comparison. Although the sample size remains 
small, the lenses acceptably reproduce the distribution of size and 
luminosity for defined by the whole a > 300 km galaxy popu¬ 

lation. SNL-0 and SNL-2 are somewhat larger and more luminous 
than average, while SNL-1 lies among the most compact of the 
high-CT galaxies. The properties of SNL-1, where leasing requires 
a « 1.2 are in fact very similar to those of bl9, where Lasker et al. 
find a = 2.2±0.6 from dynamics (unless bl9 harbours an extremely 
massive central black hole). 
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In summary, we do not find support from SNELLS 
for a strong correlation of a with compactness among the 
(j > 300 km s“ population. Clearly this conclusion is heavily de¬ 
pendent on the results for SNL-l, however, and the current small 
sample cannot rule out a correlation having large intrinsic scatter. 


6 SUMMARY 

We have presented first results from the SNELLS survey for nearby 
strong-lensing ellipticals, demonstrating the success of our infra¬ 
red lEU search method with the discovery of two new 2 < 0.05 
lenses, as well as recovering the only previously-known case, 
ES0325-G004. The survey yield, i.e. identified lenses per new ob¬ 
served candidate, is broadly compatible with the 10 per cent rate 
estimated from the Ha luminosity function of Sobral et al. (2013) 
(but note that one of the new sources, like that in ES0325-G004, is 
an [O III] emitter, and hence not included in our calculations). Com¬ 
pared to distant lens samples, SNELLS selects systems in which the 
Einstein radius is a smaller fraction of the stellar effective radius, 
so that the lensing mass is more reliably dominated by the stellar 
component. 

Our lens galaxies have high velocity dispersions 
(ct > 300 kms“^) and a-enhanced abundances ([Mg/Fe] > 0.3), 
and hence are characteristic of the systems where previous works 
have inferred heavy IMFs though dynamical (Cappellari et al. 
2013), lensing (Treu et al. 2010) or spectroscopic (Conroy & van 
Dokkum 2012b) analyses. Notwithstanding these properties, the 
SNELLS galaxies have lensing-derived stellar mass-to-light ratios 
that are individually or collectively compatible with a Milky-Way¬ 
like Kroupa (2001) IMF. A heavier-than-Salpeter IMF (without 
intrinsic scatter) is disfavoured (3.5cr significance) if we assume 
the stellar populations to be old ((a) = L10±0.08±0.10), but 
cannot be excluded if we fit the ages from the currently-available 
spectra ({a) = L23±0.15 ±0.10). Very heavy IMFs with a > 2 are 
firmly rejected using either approach. The result for ES0325-G004 
(SNL-0) is consistent with our previous estimates (Smith & Lucey 
2013), despite several changes in data and methods. The mean IMF 
result obtained from SNELLS is mildly inconsistent (at the 2-3 (t 
level) with the IMF-versus-cr trend reported by Cappellari et al. 
(2013), after accounting for intrinsic scatter and our field-of-view 
selection bias. There is a larger discrepancy with respect to the 
Treu et al. (2010) analysis of the SLACS lens sample, part of which 
can be attributed to the different treatment of projected dark-matter 
contributions. The SNELLS results are also apparently in conflict 
with the bottom-heavy IMFs inferred from spectroscopic fitting by 
Conroy & van Dokkum (2012b) at similarly high cr and Mg/Fe. 

The largest source of uncertainty in each individual measure¬ 
ment of a is in the treatment of Tref, while the extent to which we 
can infer the distribution of a in the population of giant early-type 
galaxies is limited by the small sample size. In future work with 
SNELLS, we plan to observe more lens candidates, with the goal 
of doubling the sample of low-redshift lensing systems available for 
analysis. We are also pursuing homogeneous high-signal-to-noise 
optical/IR spectroscopy of the lens galaxies. The optical spectra 
will yield more reliable estimates for the stellar age and metallicity 
of the lenses, and hence Tref, obviating the need to assume purely 
old stellar populations. The red and IR spectra will cover a range of 
IMF sensitive features (including the Nal 8200-A doublet and Fe H 
Wing-Ford band) providing a stringent test of spectroscopic mod¬ 
els in cases where the low-mass stellar content is constrained by the 
lensing M/L. Eventually, we will to use these systems to compare 
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Figure 8. Velocity dispersion, tr, effective radius Reff and total luminosity 
Lj for the three SNELLS lenses, compared to the full 6dF galaxy sample. 
Dark points show all 6dF galaxies with cr > 300 km s “ ^. The three targets 
span most of the range in size present among high-cr galaxies. For compar¬ 
ison, we show also the compact galaxy discussed by Lasker et al. (2013), 
which has structural properties similar to those of SNL-1. 

multiple methods of IMF estimation (lensing, dynamical, spectro¬ 
scopic) for a common sample of nearby galaxies, where high qual¬ 
ity data can be readily obtained. 

At present, SINFONI is the only IFU operating with suffi¬ 
ciently wide contiguous field-of-view to perform a study of this 
kind efficiently in the IR. However, the success of our search strat¬ 
egy suggests it should be worthwhile to search for lensed back¬ 
ground sources in large optical IFU datasets, e.g. from the SAMI 
(Allen et al. 2015) and MaNGA (Bundy et al. 2015) surveys. 
MaNGA is the most promising of these options, since its wider 
spectral coverage (to 1 /rm) will detect background Ha sources out 
to 2 « 0.5 and [O II] to 2 « 1.5. In the longer term, the HARMONI 
first-light instrument (Thatte et al. 2014) on the 39m European 
Extremely Large Telescope, with 5xl0arcsec^ field-of-view IFU, 
should enable such deep observations that strongly-lensed back¬ 
ground line emitters could be detected for any chosen massive el¬ 
liptical. 
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